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ABSTRACT 

We use a combination of three large A^-body simulations to investigate the depen- 
dence of dark matter halo concentrations on halo mass and redshift in the Wilkinson 
Microwave Anisotropy Probe year 5 (WMAP5) cosmology. The median relation be- 
tween concentration and mass is adequately described by a power-law for halo masses 
in the range 10^^ — 10^^ /i~^Mq and redshifts z < 2, regardless of whether the halo 
density profiles are fit using NFW or Einasto profiles. Compared with recent analyses 
of the Millennium Simulation, which uses a value of ag that is higher than allowed by 
WMAP5, z = halo concentrations are reduced by factors ranging from 23 per cent 
at 10^^ /i^^Mq to 16 per cent at 10^'' h~^MQ. The predicted concentrations are much 
lower than inferred from X-ray observations of groups and clusters. 

Key words: galaxies: clusters: general - galaxies: haloes - cosmological parameters 
- dark matter - cosmology: theory - methods: A'^-body simulations 



1 INTRODUCTION 

Investigations into dark matter (hereafter DM) clustering 
and dynamics have progressed greatly thanks to large- 
scale cosmological simulations run with the aid of super- 
computers. In particular, the A'^-body technique has al- 
lowed us to follow the detailed hierarchical build-up of viri- 
alised DM structures, resulting in near spherical haloes that 
are well described by the Navarro, Frenk & White pro- 
file (jNavarro et al.lll997l . hereafter NFW). The NFW density 
profile is described by just two parameters, the concentra- 
tion, c, and the total mass, M, of the halo. Simulations have 
shown that these two parameters are correlated, with the 
average concentration of a h alo being a weakly decreasing 
funct ion of mass (e.g NFW: iBullock et al.ll200ll: lEke et al'l 
I2OOII . hereafter ENS: IShaw et al.ll2006l : iMaccio et al.ll2007h . 
Observations of galaxy groups and clusters using X-ray and 
gravitational lensing data are being used to test this hy- 
pothesis and generally confirm t he anti-correlation betw een 
c and M (e.g. lBuote et aL I l2007l : ISchmidt fc Allenll2007l ) , al- 
though the concentrations inferred from strong gravitational 
lensing may exceed those predicted by the simulat ions (e.g. 
iHennawi et allbOOTl : iBroadhurst fc Barkanall2008l ). 

The best statistical constraints on the halo concen- 
tration distri bution currently com e from the Millennium 
Simulation ISpringel et al.ll2005^ . hereafter MS, in which 
2160^ DM particles were allowed to interact gravitationally 
in a cosmological box of length 500/i~^Mpc. The resulting 
cjM) relation is well described by a power law, c oc M"'''^ 
iNeto et al.ll2007l . hereafter N07). The MS used cosmological 



parameters from the Wilkinson Microwave Anisotropy Probe 
year 1 (WMAPl) data release (jSpergel et al.1 bOOsI ) com- 
bined with constr aints from the Two Degree Field Galaxy 
Redshift Survey (jPercival et al.ll200ll ). henceforth referred 
to as the WMAPl cosmology. 

In this letter, we use a set of A''-body simulations to 
quantify the relation between concentratio n and mass in 
the m ore recent WMAP year 5 cosmology (jKomatsu et al.l 
|2008| ). henceforth WMAP5, in which the most signifi- 
cant change from WMAPl is a downward shift in erg 
by around 13 per cent. Because the c(M) relation is 
very sensitive to the primordial power spectrum (e.g . 
ENS; 'Ala m. Bullock, fc Weinberg |2002| : iDolag et al.1 12004 
Ikuhlcn c t al.ll2005l )." it is not clear whether the MS results 
accurately describe our Universe. Indeed, we will show that 
the WMAP5 cosmology results in a significant decrease in 
both the slope and the normalisation of the c(M) relation. 
While this reduction in halo concentrations may make it 
easier for models of galaxies to match observations (e.g. 
Ivan den Bosch. Mo. fc Yand l2003l : ICnedin et alll2007l ). we 
will show that it results in strong disagreement between 
simulations and X-ray observations of groups and clusters 
of galaxies. 



2 SIMULATIONS 

We analyse a set of thre e A-body simulations, run using 
Gadget2 (|Springelll200d ). Each simulation contains 512^ 
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DM particles but with a progressively larger comoving box 
size: 25, 100, and 400 /i"^Mpc for runs L025, LlOO and L4OO, 
respectively. By combining the results of the three simula- 
tions, we cover four orders of magnitude in halo mass, a 
range that exceeds that of N07, although our total num- 
ber of haloes is smaller than in the MS. Our highest resolu- 
tion simulation (L025) uses a Plummer-equivalent comov- 
ing softening of 2kpc/h, with a maximum proper value of 
0.5kpc/h, reached at 2 = 3; LlOO and L4OO have values 
4 and 16 times larger respectively. The particle masses are 
8.34x10*^ h-^MQ, 5.33x10® h'^MQ and 3.41x10^° H'^Mq, 
for L025, LlOO and L4OO, respectively. 

Glass-like cosmological initial conditions were gener- 
ated at redshift z — 127 using the Zeldovich approximation 
and a transfer funct i on ge nerated using CMBFAST (v. 4.1, 
ISeliak fc Zaldarriagal Il996h . We use the WMAP5 (CMB 
only) cosmology, with [Qm, ^b, ^a, h, erg, rig] given by [0.258, 
0.0441, 0.742, 0.719, 0.796, 0.963]. For comparison, the MS 
used the WMAPl cosmology, [0.25, 0.045, 0.75, 0.73, 0.9, 
1.0]; for which the value of erg is about 13 per cent larger. 

2.1 Halo definitions and density profiles 

Halo virial masses and radii are determined using a 
spheric al overdensity rout ine within the SUBFIND algo- 
rithm (|Springel et al.l I2OOII ) centre d on the main su bhalo 
of Friends-of-Friends (FOF) haloes (|Davis et al.lll985l '). 

We perform all calculations for three different halo def- 
initions, which all take the halo centre to be the location of 
the most bound particle in the FOF group. According to the 
first definition, which is motivated by the spherical top-hat 
collapse model, a halo consists of all matter within the ra- 
dius Tvir for which the mean internal density is A times the 
critical density pcrit = 3H^ /SivG, where A depends on both 
cosmo logy and redshift and is given by iBrvan fc NormanI 
(jl998f l. Using the second (third) definition, a halo consists of 
all matter within the radius r2oo (''mean) for which the mean 
internal density is 200 times the critical (mean background) 
density. Note that NFW adopted the second definition. We 
will use Mvir, M200, and Mmcan to denote the corresponding 
halo masses. 

Two different density profiles (NFW and Einasto) are 
fit to well-resolved haloes using a procedure which closely 
follows N07. For each halo with at least 10* particles within 
rvir, all particles in the range —1.25 ^ logjQ(r/rvir) ^ are 
binned radially in equally spaced logarithmic bins of size 
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-0.078. Density profiles are then fit to these bins 



by minimising the difference between the logarithmic den- 
sities of the model and the data, assuming equal weighting. 
H aloes are only used if the convergence radius, as proposed 
bv lPower et"aLl lj2003l ). is within the minimum fit radius of 
0.05 Tvir. The convergence radius is defined such that the 
two-body dynamical relaxation timescale of the particles in- 
ternal to this point is similar to the age of the universe. The 
best-fit density profiles are used to define the halo radii, with 
the corresponding masses found by integrating the profiles 
to these radii. 

We consider two samples. Our default sample contains 
all haloes that satisfy our resolution criteria while our 're- 
laxed' sample retains only those haloes for which the sep- 
aration between the most bound particle and the centre of 
mass of the FOF halo is smaller than 0.07rvir. Note that 
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Figure 1. The relation between the Einasto power-law index a 
and the peak height parameter u for stacks of haloes of our default 
sample. Data from redshifts z = 0,1 and 2 was combined. Data 
points correspond to medians and error bars indicate quartiles. 
The solid curve is the best-fit quadratic function, whereas the 
dashed curve is the result from G07. The parameters for both fits 
are given in the legend together with the errors (68.2 per cent 
confidence level) obtained by bootstrap resampling the haloes. 



N07 found that this simple criterion resulted in the removal 
of the vast majority of unrelaxed haloes and as such we do 
not use their additional criteria. At z — our default and 
relaxed samples contain 1269 and 561 haloes in total. 
The NFW density profile is given by 

5c 



p{r) = Pcrit 



(r/rs)(l + r-/r,)2 



(1) 



where (5c is a characteristic density contrast and rs is a scale 
radius. The concentration is defined as C200 = '"200 /''s- NFW 
profiles were fit using the two parameters rs and 5c- 

The Einasto profile is a rolling power-law first intro- 
duced to describe the distribution of old stars in the Milky 
Way. It takes the form 



din p 
din r 



= -2 



r 
r-2 



(2) 



where r_2 is the radius at which the logarithmic slope of the 
density is isothermal (i.e. —2), analogous to in the NFW 
profile. As a result the concentration, C200 = '"2oo/r-2, is 
close to the NFW value. Integrating Eq. [2] gives 

ln(p/p_2) = --[(r/r_2)"-l] , (3) 
01 

where /9_2 is the density at r_2. iGao et all l)2007l . hereafter 
G07) have shown that there exists a simple relation between 
a and v, the dimensionless 'peak-height' parameteiQ, inde- 
pendent of redshift. We check this result for our simulations 
by performing a three-parameter fit to profiles averaged over 
ten haloes (to remove the effects of substructure). As shown 
in Fig. [T] our results are in agreement with G07; we there- 
fore adopt their a(v) relation to reduce the model to two 
free parameters (p_2 and r-2)- 



^ The peak height is defined to be ^ = <5(.rit/''"(A^viri where 
"^crit = 1.686 is the linear density collapse threshold at the present 
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Figure 2. Concentration-mass relations for z = (top), 1 (mid- 
dle) and 2 (bottom) using NFW density profiles. Data points 
correspond to median values and error bars to quartiles. Only 
bins containing at least 5 haloes are shown. The black, solid lines 
show the best-fit power-law relation. The errors on the best-fit 
parameters, given in the legend, are determined by bootstrap re- 
sampling the haloes and correspond to 68.2 per cent confidence 
limits. The pink dashed line in the top panel shows the best-fit 
power-law relation to all haloes from N07. The other curves rep- 
resent the prescriptions discussed in scction r3.il The brown solid 
lines correspond tothe special case where we set C = in Eq. 4, 
and fit a power-law to the data at each individual rcdshift. 



3 RESULTS 

Fig. [2] shows the median C2oo(A^20o, -z) relation for z = 0, 1 
and 2 (top, middle and bottom panels, respectively) using 
NFW density profiles. Tfie data, which spans four orders of 
magnitude in mass and a redshift range z — — 2 are well 
described by a function of the form 

c = A{M/Mpi.otf{l + zf. (4) 

The black solid line indicates the best-fit relation of this 
form to the data. The parameter values and errors are given 
in Table[T]for three different halo definitions, for both NFW 
and Einasto profiles, for the full and relaxed samples, and 
for both 2 = (in which case we set C = 0) and z — 0-2. 

day and cr(Mvir, z) is the linear rms density fluctuation at redshift 
z in a sphere containing a mass M^ir- 



All fits use Mpivot — 2x 10^^ h~^MQ, the median halo mass, 
in order to minimise the covariance between A and B. 

Concentration is a decreasing function of both mass 
and redshift, regardless of the model density profile and 
halo definition that is used. We note that at high redshift, 
2 = 2, the mass de pendency is significantly reduced (see 
also lZhao et al.|[2003l ). Compared with NFW concentrations, 
Einasto concentrations decrease significantly more rapidly 
with both mass and redshift. At z = the two concentra- 
tions agree for M ~ W^"^ H'^Mq. Concentrations are very 
sensitive to the halo definition. In particular, at 2 = val- 
ues for Cnioan are nearly twice as large as those for C2oo- As 
expected, the difference becomes smaller at high redshift 
because the critical density will approach the mean density 
as the redshift becomes high enough for matter to dominate 
over vacuum energy. At 2 = the median concentrations are 
typically about 10 per cent greater for the relaxed sample 
than for the default sample. 

Using C = and Mpivot = IO^/i^^Mq, N07 found 
as the best fit for aU NFW haloes in their z = sample 
(^200,-8200) ~ (4.67,-0.11). Fitting the same function to 
our 2; = sample gives (^4200,^200) ~ (3.93, —0.097) which 
yields concentrations that are lower by a factor ranging from 
23 per cent at lO" H'^Mq to 16 per cent at lO" /i^^Mq. 
This difference can be attributed to the decrease in erg. If erg 
is higher then haloes of a given mass form earlier. This in- 
creases the concentration because c refiects the background 
density of the universe at the time when the halo forms 
(NFW). 

The scatter in the concentration about the median c(M) 
relation ha s been shown to be lognormal for relaxed haloes 
IJindbOOd ). with a slight decrease in the scatter as a func- 
tion of mass (N07). The inclusion of unrelaxed haloes adds 
a tail towards low concentrations. Fig. [3] shows histograms 
of the 2 = concentrations using both NFW (black) and 
Einasto (red) profiles for both the default (solid) and re- 
laxed (dashed) samples. Either density profile yields a dis- 
tribution which agrees qualitatively with that found by N07 
for the NFW model. Fitting lognormal functions to the 
probability density functions yields (T(logjQ C200) = 0.15 and 
0.17 for the NFW and Einasto density profiles, respectively. 
For the relaxed sample the scatter is significantly smaller: 
cr(logj^Q C200) = 0.11 and 0.12. 

Fig ure [4] compares our pred ictions with the observa - 
tions of iBuote etall (|2007l ') and ISchmidt fc AUenI l|2007h . 
who measured NFW Cvir concentrations from X-ray obser- 
vations of some 70 groups and clusters with a median red- 
shift of 2 = 0.1. The observationally inferred concentrations 
are significantly greater than the predicted medians and the 
discrepancy increases with decreasing mass. Using only re- 
laxed haloes (dashed lines) does not help much because it 
removes only the low-concentration tail. Observations and 
theory could be brought back into agreement if the pre- 
dicted concentrations are too low, for example because gas 
physics has not been taken into account. Alternatively, the 



^ As a check of our analysis procedure, we ran and analysed a 
miniature version of the MS with identical mass and force reso- 
lution (we used a box of 50/i~^Mpc on a side containing 216^ 
particles) and were able to reproduce N07's result to within the 
errors. 
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Table 1. Best-fit parameters for the median c{M,z) relation (Eq. |4]l using Mpjvot = 2 X lO^'^ /i~^Mq for three different halo definitions, 
two different density profiles, redshift 2 = and z = 0—2 and for both the full (F) and relaxed (R) halo samples. The errors correspond 
to 1 CT confidence intervals and have been determined by bootstrap resampling the haloes. 
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Figure 3. Histogram of halo concentrations using the NFW 
(black) and Einasto (red) density profiles for the default (solid) 
and relaxed (dashed) samples. 



Figure 4. Comparison of observed and predicted NFW Cvir(Af) 
relations at 2 = 0.1. Data point s indicate weighte d averages 
over about 10 ob jects taken from iBuote et al. I ||2007|') (squares) 
and ISchmidt fc AUenl l|2007l ') (diamonds). The sets of solid 
(dashed) lines are our predictions for the median concentrations 
and the +lcr lognormal scatter. 



observed sample may be highly biased towards objects with 
high concentrations because these typically have higher X- 
ray luminosities, particularly for groups. 



3.1 Comparison with predictions from the 
Uterature 

A number of semi-empirical prescriptions (i.e. theoretically 
motivated fitting functions) have been proposed to predict 
the c(M, z) relation for arbitrary cosmologies. The basic 
premise of each is that the concentration of a halo reflects 
the background density of the universe at the formation time 
of the halo. The prescriptions differ mainly in the definition 



of this formation time. We will compare our results to the 
prescriptions proposed by NFW and ENS0 

The NFW prescription for c(M) has three free parame- 
ters. The first two, F and /, are used to define the collapse 
redshift as the redshift for which a fraction F of the final 
halo mass M is contained in progenitors of mass ^ fM. 
The third parameter, C, gives the ratio of 5c, which is 1-1 
related to the concentration parameter in the NFW prescrip- 
tion, and the mean density of the universe at the collapse 
redshift. The values initially suggested by NFW are [F, f, 



2 We also compared with lBu llock et al.l ll20 0l | ). using the values 
of K = 2.9 and F = 0.001 jWechsler et al.ll2006l) and found sig- 
nificant disagreement with our data at high masses (see Fig. |2]l. 
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C) = (0.5, 0.01, 3000), but we find that the modified values 
(0.1, 0.01, 200) provide a much better fit to our default sam- 
ple while the relaxed sample is best fit by (0.5, 0.01, 500) 
although the latter values are poorly constrained. Both the 
original and our modified prescriptions are compared to our 
data in Fig. [5] The modified NFW model reproduces the 
full sample of NFW halo concentrations over a wide range 
of masses and redshifts, while the original model strongly 
overestimates the concentrations at z = 2 (bottom panel). 

The ENS prescription has only one free parameter, 
Cct, which implicitly defines the collapse redshift through 
_D(zc)o"efl(Ms) = 1/Ccr, where D{z) is the linear growth fac- 
tor, aeff(Af) is the effective amplitude of the linear power 
spectrum at z = and Ms = M{r < 2.17rs). The concentra- 
tion then follows (iteratively) by equating the characteristic 
density ps = 3M/{ATvrl) to the spherical collapse top-hat 
density at the collapse epoch. ENS found that Co- ~ 28 was 
a good fit, but G07 concl uded that this did not describe their 
data very well (see also lZhao et al.ll2003l ). We performed a 
fit (equally weighting all points) and found Co- ~ 29 and 
30 as the best-fit values using NFW and Einasto density pro- 
files, respectively. While the ENS prescription is an excellent 
fit for NFW density profiles below z = 2, it underestimates 
the 2 = concentrations for M < lO^'^ h~^MQ when the 
Einasto density profile is used. 



4 CONCLUSIONS 

We have measured the concentrations of DM haloes us- 
ing three cosmological A''-body simulations assuming the 
WMAP5 cosmology. We presented power-law relations be- 
tween halo concentration and mass for z = and z = — 2 
using both NFW and Einasto density profiles for three dif- 
ferent halo definitions and we compared them to predic- 
tions from the literature. We found that halo concentrations 
are significantly lower in the WMAP5 cosmology than in 
the WMAPl cosmology, which was used in the Millennium 
Simulation. For 2 = the reduction varies from 23 per cent 
at 10" ft-^Mg to 16 per cent at 10^"* H'^Mq. 

While the decrease in the concentrations may improve 
the agreement between models and observations of galaxies, 
we found that it results in significant discrepancies with X- 
ray observations of groups and clusters of galaxies. To deter- 
mine the seriousness of this discrepancy, it will be necessary 
to carefully study possible observational selection biases as 
well as the effects of baryons on the halo concentrations. 
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